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ABSTRACT: This review surveys the observed properties of interstellar dust grains: the 
wavelength-dependent extinction of starlight, including absorption features, from UV to infrared; 
optical luminescence; infrared emission; microwave emission; optical, UV, and X-ray scattering 
by dust; and polarization of starlight and of infrared emission. The relationship between presolar 
grains in meteorites and the interstellar grain population is discussed. Candidate grain materials 
and abundance constraints are considered. A dust model consisting of amorphous silicate grains, 
graphite grains, and polycyclic aromatic hydrocarbons is compared with observed emission and 
scattering. Some issues concerning evolution of interstellar dust are discussed. 
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1 INTRODUCTION 



Dust grains play a central role in the astrophysics of the interstellar medium, from 
the thermodynamics and chemistry of the gas, to the dynamics of star formation. 
In addition, dust shapes the spectra of galaxies: radiation at short wavelengths 
is attenuated, and energy is radiated in the infrared. It is estimated that 30% or 
more of the energy emitted as starlight in the Universe is reradiated by dust in 
the infrared (Bernstein et al. 2002). Interstellar dust determines what galaxies 
look like, how the interstellar medium (ISM) in a galaxy behaves, and the very 
process of star formation that creates a visible galaxy. 

This review summarizes current knowledge of the abundance, composition, 
and sizes of interstellar dust grains, as indicated by observations of extinction, 
scattering, and emission from dust, supplemented by evidence from interstellar 
gas phase abundances and presolar grains in meteorites. The term "dust grain" is 
understood here to extend down to molecules containing tens of atoms, as there 
is no discontinuity in the physics as the particle size decreases from microns to 
Angstroms. 

The distribution of gas and dust in the Galaxy - the structure of the interstellar 
medium - is not discussed. 

The astrophysics of interstellar dust is not covered in this review. This includes 
the physical optics of small particles, charging of dust grains, heating and cooling 
of dust grains, chemistry on dust grain surfaces, forces and torques on dust grains, 
sputtering and shattering of dust grains, and alignment of dust with the magnetic 
field; an introduction to these topics can be found elsewhere (Kriigel 2002, Draine 
2003a). 

The review does not attempt to survey or compare different models proposed 
for interstellar dust. One dust model is introduced to illustrate modeling of 
absorption, scattering, and infrared emission from interstellar dust. 

It is not possible to cite all of the important papers in this area; only a few 
articles are cited in connection with each topic. The reader is encouraged to 
also consult prior reviews by Savage & Mathis (1979), Mathis (1990, 1993, 2000), 
Draine (1995), Witt (2000a), Witt (2000b), Voshchinnikov (2002), and the book 
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by Whittet (2003). Dorschner & Henning (1995) discuss the metamorphosis of 
interstellar dust. 

2 EXTINCTION 

The existence of interstellar dust was first inferred from obscuration, or "extinc- 
tion", of starlight (Trumpler 1930). Much of our knowledge of interstellar dust 
continues to be based on studies of the wavelength-dependence of this attenua- 
tion, often referred to as "reddening" because of the tendency for the extinction to 
be greater in the blue than in the red. The wavelength-dependence strongly con- 
strains the grain size distribution, and spectral features (§3) reveal the chemical 
composition. 

The extinction is most reliably determined using the "pair method" - com- 
paring spectrophotometry of two stars of the same spectral class; if one star has 
negligible foreground dust while the second star is heavily reddened, compari- 
son of the two spectra, together with the assumption that the dust extinction 
goes to zero at very long wavelength, allows one to determine the extinction 
A\ = 2.51og 10 (i ? x/i ; A) as a function of wavelength A, where F\ is the observed 
flux and F® is the flux in the absence of extinction. The "pair method" has been 
used to measure extinction curves for many sightlines, in many cases over a range 
of wavelengths extending from the near-infrared to the vacuum UV. 

2.1 Milky Way Dust 

2.1.1 OPTICAL-UV EXTINCTION CURVES 

The dimensionless quantity Ry = Ay/(As — Ay) is a common measure of the 
slope of the extinction curve in the optical region. Very large grains would pro- 
duce gray extinction with Ry —>■ oo. Rayleigh scattering (A\ oc A -4 ) would 
produce very steep extinction with Ry ~ 1.2. Ry is known to vary from one 
sightline to another, from values as low as 2.1 (toward HD 210121) (Welty & 
Fowler 1992) to values as large as 5.6-5.8 (toward HD 36982) (Cardelli et al. 
1989, Fitzpatrick 1999). Cardelli et al. (1989) showed that normalized extinction 
curves A\/Aj (using the / band extinction to normalize) could be approximated 
by a seven-parameter function of wavelength A: 

A x /Aj « /(A; Ry,C 1 ,C 2 , C 3 ,C 4 ,A ,7) • (1) 

At wavelengths A > 3030A, the function /(A) depends only on A and the single 
parameter Ry. The parameters C3, Ao, and 7 determine the strength and shape 
of the 2175A "bump", and the coefficients C\, C2, and C4 determine the slope 
and curvature of the continuous extinction at A < 3030A. 

Cardelli et al. (1989) show that if the single quantity Ry is known, it is possi- 
ble to estimate the values of the other six parameters (Ci_4, Ao, 7) so that the 
optical- UV extinction can be approximated by a one-parameter family of curves. 
However, if the UV extinction has been measured, an improved fit to the obser- 
vations can be obtained by fitting Ci_4, Ao and 7 to the measured extinction. 

Empirical extinction curves in diffuse clouds show relatively little variation in 
the infrared, and for wavelengths 0.7/xm £ A ^ 8/xm the function (1) appears 
to be approximately "universal" (i.e., independent of Ry) in diffuse clouds [in 
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Figure 1: Extinction curves from prescription of Fitzpatrick (1999), with diffuse interstellar 
bands (DIBs) added as described in §3.3. The DIBs are barely visible on this plot. 



dense clouds grains acquire "ice" mantles (§3.5) which alter the extinction]. Fitz- 
patrick (1999) pays careful attention to the effects of finite-width photometric 
bandpasses, and gives modified formulae which appear to improve the overall fit 
to observations of ice-free dust. The Cardelli et al. (1989) and Fitzpatrick (1999) 
fits for Ry = 3.1 are compared in §8 (Table 4). 

Figure 1 shows extinction curves calculated using the F99 parametrization for 
Ry = 2.1, 2.5, 3.1, 4.0, and 5.5. The coefficients in the Cardelli et al. (1989) or 
Fitzpatrick (1999) fitting formulae can be adjusted to improve the fit to specific 
sightlines; such a fit is shown for the extreme case of HD 210121, showing that 
the UV extinction can differ significantly from the average behavior for the same 
value of Ry. 

"Pair method" determinations of the reddening law for many sightlines indicate 
that Ry ~ 3.1 for the "average" extinction law for diffuse regions in the local 
Milky Way (Savage & Mathis 1979, Cardelli et al. 1989). Sightlines intersecting 
clouds with larger extinction per cloud tend to have larger values of Ry] the 
larger Ry values may indicate grain growth by accretion and coagulation. 

Another approach to determining the wavelength-dependent extinction is to use 
star counts or galaxy counts as a function of apparent brightness. Szomoru & 
Guhathakurta (1999) use UBVRI photometry of Galactic field stars to determine 
the extinction law; for the four high-latitude clouds they studied, they conclude 
that Ry ;$ 2, well below the value Ry « 3.1 which is widely considered to be 
"average" , and comparable to the most extreme values of Ry ever found in studies 
of individual stars (Ry = 2.1 for HD 210121). Further study is needed to reconcile 
this apparent conflict. 

Udalski (2002) uses V and / photometry of "red clump giants" to study the 
reddening law. Toward the LMC he infers Ry ~ 3.1, but toward Galactic bulge 
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regions he finds Ry « 1.8 — 2.5. 

Sloan Digital Sky Survey (SDSS) photometry (Lupton et al. 2003) can be 
used to study the reddening toward stars at the "blue tip" of the main sequence 
(Finkbeiner et al. 2003b). The observed reddening appears to be consistent with 
Ry ss 3.1. Galaxy surface brightnesses can also be used; SDSS galaxy photometry 
is consistent with Ry = 3.1 (Schlegel et al. 2003). 

2.1.2 EXTINCTION PER H 

Using H Lyman-a and absorption lines of H2 to determine the total H column 
density Ah, Bohlin et al. (1978) found 

Nh/(A b - Ay) = 5.8 x 10 21 cm-Vag- 1 , (2) 

N H /A V « 1.87 x 10 21 cm-Vag- 1 for Ry = 3.1 , (3) 

to be representative of dust in diffuse regions. For Ry = 3.1, the F99 reddening 
fit gives Aj c /Ay = 0.554 for Cousins I band (A = 0.802 /mi), thus 

A Ic /N H « 2.96 x 10" 22 magcm 2 for Ry = 3.1 . (4) 

If the extinction curve is really universal for A ^ 0.7 /xm we might expect eq. (4) 
to apply independent of Ry. 

Rachford et al. (2002) report H and H2 column densities for 16 sightlines 
through diffuse and "translucent" clouds. Figure 2 shows A[ c /Nn versus Ry 1 for 
the 14 sightlines with Ah known to better than a factor 1.5 . The Ic band extinc- 
tion Aj c has been estimated from Ab — Ay and Ry using the F99 i2y-dependent 
reddening law. 

Figure 2 appears to show that Ai c /Nu is not universal for local Milky Way 
dust: the Ic band extinction per H deviates from the ratio (4) when Ry deviates 
from 3.1, in the sense that Aj c /Nn increases when Ry is larger (larger grains). 
We note that the trend in Aj c /Nn could be in part the result of errors in de- 
termination of Ry, (i.e., errors in estimating Ay, since E(B — V) is probably 
accurate to better than 10% for these stars), but the data appear to show a real 
trend. A least-squares fit, constrained to pass through (4), gives 

A Ic /N H » [2.96 - 3.55(3.1/i?y - 1)] x 10~ 22 magcm 2 . (5) 

Figure 3 shows A\/Nk, the extinction per H, estimated from the F99 extinction 
fits and eq. (5). If the F99 extinction fits and eq. (5) are both correct, the 
extinction per H increases for A ^ 0.3 fim when Ry increases. Presumably this 
is due to increased scattering by the grains as the grain size increases. 

2.1.3 INFRARED EXTINCTION 

Between ~ 0.9 fim and ~ 5^m the continuous extinction curve can be approx- 
imated by a power-law, A\ oc \~@, with [5 ~ 1.61 (Rieke & Lebofsky 1985), 
1.70 (Whittet 1988), 1.75 (Draine 1989b), ~1.8 (Martin & Whittet 1990), ~1.8 
(Whittet et al. 1993), or 1.70 (Bertoldi et al. 1999, Rosenthal et al. 2000). 

The extinction in the 5-8 /Ltm region is controversial (see Figure 4) . Draine 
(1989b) concluded that the observational evidence was consistent with A\ oc 
A" 1 - 75 out to A ?3 6 /im. Using ISO observations of H recombination lines, Lutz 
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Figure 2: Ai/Nu for 14 sightlines through translucent clouds (Rachford et al. 2002), as a 
function of 1/Ry, with Ry determined by IR photometry (circles, Rv ±0.2); polarization A max 
(triangles, Rv ± 10%); or UV extinction curve (diamonds, Rv ± 10%). Vertical error bars show 
l-o" uncertainty due to errors in A?h; tilted error bars show effects of errors in Rv- Least-squares 
fit [eq. 5)] is shown. 



et al. (1996) find that the extinction toward Sgr A* does not decline with increas- 
ing A in the 4-8 am region. However, studies using H2 rovibrational lines in the 
Orion molecular cloud, (Bertoldi et al. 1999, Rosenthal et al. 2000) find that the 
extinction continues to decline with increasing A to a minimum at ~6.5 /iin. The 
wavelength dependence of the continuous dust extinction in the Orion molecular 
cloud could conceivably differ from that on the sightline to Sgr A* ; alternatively, 
perhaps the H recombination line intensity ratios differ from the "case B" re- 
combination assumed by Lutz et al. It is interesting to see in Figure 4 that the 
Weingartner & Draine (2001a) grain model for Ry = 5.5 is in rough agreement 
with the Lutz et al. results toward the galactic center. Further observational 
study of the 5-8 /xm extinction is needed. 

In addition to this continuous extinction, there are strong absorption features 
at 9.7/um and 18 urn due to silicates (§3.2), a hydrocarbon feature at 3.4/xm 
(§3.4), and PAH absorption at 6.2 am (§3.6). The extinction in molecular clouds 
also shows additional bands due to ice mantles (§3.5). 

2.2 Dust in Other Galaxies 

Dust in the Milky Way varies from one sightline to another. It is obviously 
of great interest to study the dust in other galaxies, both in order to correct 
observations for the extinction by the dust, but also to inform our understanding 
of dust astrophysics. 

Our knowledge of the extinction law for dust in other galaxies is greatest for 
the Large and Small Magellanic Clouds (LMC and SMC), where we can study the 




extinction for individual stars, but some progress has been made in determining 
the wavelength-dependence of extinction in more distant galaxies. 



2.2.1 MAGELLANIC CLOUDS 

Not surprisingly, there are regional variations within the LMC and SMC. Stars 
more than ~500 pc away from the 30 Doradus region of the LMC have Ry ~ 3.2 
and an extinction law that appears to be quite similar to the Milky Way diffuse 
cloud extinction law (Fitzpatrick 1986, Misselt et al. 1999). The reddening per 
H atom E(B - V)/N n « 4.5 x 10~ 23 cm 2 /H (Koorneef 1982, Fitzpatrick 1985), 
approximately 40% of the Milky Way value. This is approximately consistent with 
the LMC metallicity: LMC H II regions have Ne/Hw (3.5 ± 1.4) x 10" 5 (Kurt & 
Dufour 1998), ~30% of the solar Ne/H (= 1.2 x 10~ 4 ) (Grevesse k Sauval 1998). 
Stars in or near 30 Dor have extinction curves which have a weakened 2 175 A 
feature (Fitzpatrick 1986, Misselt et al. 1999). 

In the SMC, stars in the "bar" region have extinction curves which appear to 
lack the 2175A feature. The reddening per H is only E{B - V)/Nn = 2.2 x 
10~ 23 cm 2 /H, ~13% of the Milky Way value, again consistent with the SMC 
metallicity: SMC H II regions have Ne/H^ (1.6 ±0.2) x 10" 5 (Dufour 1984, Kurt 
& Dufour 1998), or 14% of the solar Ne/H, and the gas-phase C/H is only 6% of 
solar. 

Weingartner & Draine (2001a) are able to reproduce the observed extinction 
laws in the LMC and SMC using appropriate mixtures of carbonaceous and sil- 
icate grains. The sightlines in the SMC bar which lack the 2175A extinction 
feature can be reproduced by models which lack carbonaceous grains with radii 
a 0.02 /im. 
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Figure 4: Infrared extinction, relative to extinction at K (2.2 /xm), inferred by Rosenthal et al. 
(2000) for dust in the OMC-1 molecular cloud, and by Lutz et al. (1996) for dust toward Sgr 
A*, and as calculated for the Weingartner & Draine (2001a) grain models for Rv = 3.1 (60ppm 
C in PAH) and Rv = 5.5 (30ppm C in PAH). The 6.2 /jm absorption feature in the Weingartner 
& Draine (2001a) grain models is due to PAHs. 

2.2.2 M31 

Bianchi et al. (1996) present UV extinction curves toward selected bright stars 
in M31, finding an extinction law similar to the average Milky Way extinction 
curve, but with the 2175A feature possibly somewhat weaker than in the Milky 
Way. 

2.2.3 OTHER GALAXIES 

There are several different approaches to determining the extinction curve for 
dust in galaxies where individual stars cannot be resolved. 

The dust extinction law can be determined when a foreground galaxy overlaps 
a background galaxy. Berlind et al. (1997) find the extinction curve for dust 
in the spiral arm of NGC 2207 to resemble the Milky Way extinction law for 
Ry = 5.0; the interarm dust appears to be even grayer. However, the flatness 
of the extinction curve could be due in part to unresolved optically-thick dust 
patches in the foreground galaxy. 

Keel & White (2001a) find the extinction law in the spiral galaxy AM 1316-241 
to be close to the Milky Way mean extinction, with Ry = 3.4 ±0.2; for the spiral 
galaxy AM 0500-620 the dust extinction law appears to be somewhat steeper 
than the Milky Way average, with R v sa 2.5 ± 0.4. Keel & White (2001b) find 
the extinction law in the Sc galaxy NGC 3314A to be close to the Milky Way 
mean, with Ry = 3.5 ±0.2 over galactocentric radii 1.6-3.8 kpc, with no evidence 
for any radial trend in reddening law. 
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From observations of the overall emission spectrum Calzetti (2001) infers the 
internal extinction by dust in starburst galaxies. The inferred reddening law is 
relatively gray compared to the Galactic reddening law, and shows no evidence 
of a 2175A feature. It is not clear to what extent the flatness of the apparent 
reddening law may be due to the effects of radiative transfer in optically-thick 
distributions of stars and dust. 

Searches have been made for the 2175A extinction feature in QSO spectra 
(McKee & Petrosian 1974). Malhotra (1997) reports a statistical detection in a 
sample of 92 QSOs, but Pitman et al. (2000) argue that the detection was not 
statistically significant. 

Gravitationally-lensed QSOs with multiple images can be used to determine 
the extinction law if the intrinsic spectrum is assumed to be time-independent. 
Falco et al. (1999) determine extinction laws in galaxies out to z = 1.01, with Ry 
estimates from 1.5 to 7.2. Toft et al. (2000) find the extinction curve in a galaxy 
at z = 0.44 to have 1.3 < R v < 2.0. Motta et al. (2002) find that the dust in a 
lensing galaxy at z = 0.83 is consistent with a standard Milky Way extinction 
law for Ry = 2.1 ± 0.9, including an extinction bump at 2175A. 



3 SPECTROSCOPY OF DUST 



The composition of interstellar dust remains controversial. While meteorites (§4) 
provide us with genuine specimens of interstellar grains for examination, these 
are subject to severe selection effects, and cannot be considered representative of 
interstellar grains. Our only direct information on the composition of interstellar 
dust comes from spectral features in extinction, scattering, or emission. 

3.1 2175 A Feature: Aromatic C? 

By far the strongest spectral feature is the broad "bump" in the extinction curve 
centered at ~2175A (see Figure 1). 

Fitzpatrick $z Massa (1986) showed that the observed 2175A feature can be 
accurately fit with a Drude profile 



AC ext (A) = 7T7 , , (6) 



(A/A - Ao/A) 2 + 7 2 

peaking at Ao = 2175A (or Ag 1 = 4.60 /im™ 1 ) and with a broadening parameter 
7 = 0.216 (corresponding to FWHM 7 A = 469A, or 7A0 1 = 0.992 ^m" 1 ). The 
strength of the absorption is such that on an average diffuse cloud sightline in 
the Milky Way, the 2175A profile corresponds to an oscillator strength per H 
nucleon nxfx/nu ~ 9.3 x 10~ 6 (Draine 1989a). Since we expect the oscillator 
strength per molecule fx ~ 0.5, it can be concluded that the carrier molecule 
must contain one or more elements from {C,0,Mg,Si,Fe}, as these are the only 
condensible elements with sufficiently high abundances kx/jih ~4x 10 -5 in the 
ISM. 

Fitzpatrick & Massa (1986) show that the width of the feature varies consider- 
ably from one sightline to another. The FWHM has an average value 0.992 fira^ 1 , 
but the distribution has a ±2a variation of ±12%. This variability in FWHM is 
in contrast to the fact that the central wavelength Ao is nearly invariant from one 
sightline to another: although variations are detected, Ao shows ±2a variations 
of only ±0.46%. 
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Immediately following discovery of the 2175A feature by Stecher (1965), Steelier 
& Donn (1965) pointed out that small graphite particles would produce absorp- 
tion very similar to the observed feature. Although non-carbonaceous carriers 
have been proposed (e.g., OH~ on small silicate grains: Steel & Duley 1987), 
it now seems likely that some form of graphitic carbon is responsible for the 
observed feature. 

The C atoms in graphite are bound to one another in hexagonal sheets. Each C 
atom has 4 valence electrons (2s 2 2p 2 ); 3 of these are in trigonal sp 2 or a orbitals; 
the remaining electron is in a delocalized tt orbital, which is shared among the C- 
C bonds. The individual sheets of carbon atoms are held parallel to one another 
only by weak van der Waals forces. Absorption of a photon can excite a tt electron 
to an excited orbital (vr*); this tt — > tt* transition is responsible for the absorption 
feature peaking at ~ 2175A. Using the dielectric function of graphite (Draine & 
Lee 1984), one finds that randomly-oriented small graphite spheres would have 
an oscillator strength / = 0.16 per C atom (Draine 1989a). Thus the observed 
2175A extinction feature would require C/H=5.8 x 10~ 5 in small graphite spheres, 
compatible with interstellar abundance constraints. 

Although graphite is an attractive candidate, the graphite hypothesis does not 
appear to have a natural way to accomodate the observed variations in FWHM 
of the profile while at the same time holding Ao nearly constant: variations in 
graphite grain shape and size produce variations in FWHM, but these are ac- 
companied by changes in Ao (Draine & Malhotra 1993). 

Mennella et al. (1998) propose that the 2175A feature is due to UV-processed 
HAC particles, but do not appear to have actually measured in the lab a profile 
matching the interstellar 2175A feature. 

Wada et al. (1999) propose that the 2175A feature is due to onion-like hyper- 
fullerene carbon particles. While these particles do have an absorption peak near 
2200A, the feature appears to be broader than the interstellar profile. 

The carbon atom skeleton of polycyclic aromatic hydrocarbon (PAH) molecules 
is very similar to a portion of a graphite sheet, with similar electronic wavefunc- 
tions. It is therefore not surprising that PAH molecules generally have strong 
TT — > TT absorption in the 2000-2500A region. Thus large PAH molecules are 
candidates to be the carrier of the interstellar 2175A feature - this is a natural 
extension of the graphite hypothesis. 

The grain model of Weingartner & Draine (2001a) and Li &: Draine (2001b) 
has C/H=6.0 x 10~ 5 in PAH molecules or clusters containing from 20 to 10 5 C 
atoms - this population of PAH molecules is required to reproduce the observed 
infrared emission (see §9). Absorption profiles are not known for PAH molecules 
of the sizes characteristic of the interstellar PAHs, but one would expect a similar 
tt — > tt* oscillator strength per C as for graphite, or / « 0.16. Thus we see that 
in this grain model, the 2175A extinction feature is expected to be primarily - 
perhaps entirely - due to large PAH molecules. 

In this interpretation, the observed band profile would be due to a mixture of 
PAHs, and the observed variations in FWHM (and small variations in Ao) would 
result from differences in the PAH mix from one sightline to another. 

3.2 Silicate Features 

The infrared extinction includes a strong absorption feature peaking at ~9.7^m. 
Silicate minerals generally have strong absorption resonances near 10 /j,m due 
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to the Si-0 stretching mode, and it seems virtually certain that the interstellar 
9.7 /im feature is due to absorption by interstellar silicate material. This con- 
clusion is strengthened by the fact that a 10 /xm emission feature is observed 
in outflows from cool oxygen-rich stars (which would be expected to condense 
silicate dust) but not in the outflows from carbon-rich stars (where silicates do 
not form, because all of the oxygen is locked up in CO). There is also a broad 
feature at 18 /im which is presumed to be the O-Si-0 bending mode in silicates 
(McCarthy et al. 1980; Smith et al. 2000). 

The Trapezium (Gillett et al. 1975) emission implies r(A) with FWHM«3.45 /im, 
whereas ft Cep (Russell et al. 1975) has FWHM=2.3 / um, so the silicate properties 
are not universal. On sightlines dominated by relatively diffuse clouds (e.g., to- 
ward the B5 hypergiant Cyg OB2-12 or distant WC stars) the 9.7/xm absorption 
appears to be better fit by the narrower /j, Cep profile (Roche & Aitken 1984, 
Bowey et al. 1998). The silicates in molecular clouds, however, seem better fit by 
the broader Trapezium profile (Bowey et al. 1998). 

The strength of the 9.7 fim feature relative to Ay has been measured toward 
Cyg OB2-12 and toward WC stars. The studies in Table 1 together indicate 
Ay/Arg^ = 18.5 ±2.0, where the uncertainty includes a subjective allowance for 
systematic errors in the model-fitting. 
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The interstellar 9.7 //m feature is broad and relatively featureless, as opposed to 
absorption profiles measured in the laboratory for crystalline silicates, which show 
considerable structure which can be used to identify the precise mineral. The 
absence of substructure in the interstellar profile is believed to indicate that the 
interstellar silicates are largely amorphous rather than crystalline. Amorphous 
or glassy silicates can be produced in the laboratory by ion bombardment of 
initially crystalline material (Kraetschmer & Huffman 1979), formation in smokes 
(Day 1979), rapid quenching of a melt (Jager et al. 1994) or deposition following 
evaporation (Koike & Tsuchiyama 1992, Stephens et al. 1995, Scott & Duley 
1996). The laboratory absorption profiles are in some cases quite similar to the 
interstellar profile, supporting the view that interstellar silicates are amorphous. 

The observed strength per H nucleon of the interstellar silicate band appears to 
require that close to 100% of the solar abundance of Si, Fe, and Mg be condensed 
into amorphous grains. Mathis (1998) has argued that the observed strength can 
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only be reproduced if silicate grains are "fluffy", with >25% vacuum, as in the 
grain model of Mathis & Whiffen (1989). 

The observed 10 //m profile appears to be consistent with amorphous material 
with a composition in the olivine family, Mg2 X Fe2-2xSi04. Since Mg and Fe are 
of approximately equal abundance in the Sun, it would be reasonable to have 
x ps 0.5. Li & Draine (2001a) conclude that spectrophotometry of the 9.7/tm 
feature in extinction toward Cyg OB2-12 (Bowey et al. 1998) limits the fraction 
of the Si in crystalline silicates to be ^, 5%. Demyk et al. (1999) find that for 
two protostellar objects, at most 2% of the silicates could be crystalline. 

However, Bowey & Adamson (2002) argue that the observed 10 /tm profile 
toward Cyg OB2-12 is consistent with a mixture of amorphous silicates (40% by 
mass) and 8 crystalline pyroxenes (60% by mass) ; each of the crystalline species in 
their mix contributes ~7.5% by mass, and the fine structure in the 10 /im profile 
appears to be consistent with the observations. However, this interpretation 
raises various concerns: 

• The lab measurements were for powders in a KBr matrix; spectra in vacuo 
would be different. Bowey & Adamson (2002) state that the "KBr shift" 
is minimal, but it would be expected to be significant near these strong 
resonances. 

• A mix of crystalline grains that gives a smooth absorption profile will not 
give a smooth emission profile unless the different crystalline materials are 
heated to the same temperature. 

• The different crystalline types have distinct far- infrared modes. The smooth- 
ness of the far-infrared spectra of star-forming regions places upper limits 
on crystalline abundances. For example, the proposed mixture has ~7.5% 
of the Si atoms in diopside, but the observations of Onaka & Okada (2003) 
(see below) limit diopside to ^ 2% of interstellar Si. 

• It is not apparent that the proposed mixture would be consistent with the 
18 /tm extinction profile. 

• The polarization profile of grains near strong resonances is sensitive to the 
dielectric function. It is not apparent that the proposed mixture would be 
consistent with spectropolarimetric constraints (see below). 

Further study of these points is urgently required. In the text below it is assumed 
that interstellar silicates are predominantly amorphous, but the reader should 
keep in mind that a significant crystalline fraction may be possible. 

The 9.7 and 18 /tm silicate features can be observed in polarization, either in 
absorption or in emission. For strong resonances (such as the two silicate fea- 
tures), the polarization profile differs from the absorption profile, so that spec- 
tropolarimetry can constrain the dielectric function of the grain material (Martin 
1975). The observed intensity and polarization across the 9.7 /im feature appear 
to be approximately consistent with the "astronomical silicate" dielectric func- 
tion estimated by Draine & Lee (1984), but the observed 20 /mi polarization is 
stronger than expected for that dielectric function (Aitken et al. 1989, Smith et 
al. 2000). 

Onaka & Okada (2003) report detection of an emission feature at 65 /im which 
may be due to diopside CaMgSi206, containing 5-10% of interstellar Ca (cor- 
responding to 0.35-0.7% of interstellar Si); if confirmed, this would be the first 
evidence for crystalline silicates in the ISM. While crystalline silicates are not 
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abundant in the ISM, there is evidence from distinctive emission features for 
their presence in dust disks around main-sequence stars (see Artymowicz 2000 
for a summary), young stellar objects and evolved stars (see Waelkens et al. 2000 
for a summary), and in some solar-system comets (see Hanner 1999 for a sum- 
mary); see §5.1 below. Nevertheless, even in these systems most of the silicate 
material remains amorphous. 
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Figure 5: Diffuse interstellar bands in the 5700-6670 A region, using DIB parameters from 
Jenniskens & Desert (1994). 



3.3 Diffuse Interstellar Bands 

The observed extinction curve includes a large number of weak extinction fea- 
tures, known as the "diffuse interstellar bands" , or DIBs. The DIBs have FWHM 
£ lA, too broad to be due to molecules with ^ 5 atoms in the gas phase, and 
are therefore features of the grain population. Jenniskens & Desert (1994) list 
154 "certain" DIBs in the interval 0.38-0.868 fim, plus an additional 52 "prob- 
able" detections. Figure 5 shows an average extinction curve for 1.5/um -1 < 
A -1 < 1.75 jtim" 1 , a wavelength range which includes several strong DIBs, most 
notably the DIB at 5780A. This curve consists of the smooth extinction calcu- 
lated using the Fitzpatrick (1999) parametrization plus individual diffuse bands 
approximated by Drude profiles with central wavelengths, FWHM, and strengths 
from the tabulation of Jenniskens & Desert. 

The first DIBs were discovered over 80 years ago (Heger 1922) and their in- 
terstellar nature was established 69 years ago (Merrill 1934), yet to date not 
a single DIB has been convincingly identified! Some of the DIBs may be ab- 
sorption features produced by individual PAH molecules in the interstellar PAH 
mixture. Desert et al. (1995) report a positive correlation between the strength of 
the 2175A bump and the strength of certain DIBs. High resolution observations 
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of the 5797A and 6614A DIBs show fine structure which could be indicative of 
rotational bands (Kerr et al. 1996, 1998). 

DIBs have recently been detected toward reddened stars in the LMC and SMC 
(Ehrenfreund et al. 2002). 

If the DIBs were due to absorption in the aligned grains responsible for the po- 
larization of starlight (§6) there should be excess polarization associated with the 
DIBs (Martin & Angel 1974, Pureed & Shapiro 1977). All attempts to measure 
polarization structure associated with DIBs have yielded only upper limits (see 
Adamson & Whittet 1995) that appear to require that at least the DIBs that 
have been studied originate in non-aligned material, consistent with a molecular 
origin. 

At this time the carriers of the DIBs remain unknown, but it seems likely 
that at least some are due to large molecules/ultrasmall grains, possibly of PAH 
composition. 

3.4 3.4^m Feature: Aliphatic C-H Stretch 

A broad extinction feature at 3.4 /im is measurable on sightlines with Ay ^ 10 mag 
(Adamson et al. 1990). The feature appears to be present in diffuse atomic re- 
gions (in contrast to the "ice" features, see §3.5), and therefore is likely due to 
"refractory" grain material. The C-H stretching mode in aliphatic (chain-like) 
hydrocarbons occurs at 3.4 /mi. Since such hydrocarbons are plausible grain con- 
stituents, the observed 3.4 //m absorption is generally attributed to the aliphatic 
C-H stretch (Sandford et al. 1991). 

Unfortunately, aliphatic hydrocarbons show considerable variation in the C-H 
band strength, so it has not been possible to use the measured 3.4 //m feature 
to determine the abundance of aliphatic hydrocarbon material in the interstellar 
grain population, nor is it possible to identify the specific hydrocarbon mate- 
rial (Pendleton & Allamandola 2002). Sandford et al. (1991, 1995) suggest that 
the 3.4 iim. feature in the diffuse ISM is due mainly to short saturated aliphatic 
chains containing ~4% of the total C abundance. Greenberg et al. (1995) sug- 
gest that the 3.4 /im feature is due to carbonaceous organic residues produced by 
UV photolysis of ice mantles; Duley et al. (1998) attribute the feature to hydro- 
genated amorphous carbon (HAC) material incorporating 20-25% of the total C 
abundance; 

In grain models where the 3.4 /im absorption is assumed to arise in hydrocarbon 
mantles coating silicate grains (e.g., Li & Greenberg 1997), the 3.4 /im feature 
would be polarized when the 9.7/tm silicate feature is polarized; this allows a 
direct test of the silicate core/hydrocarbon mantle model for interstellar dust (Li 
& Greenberg 2002). 

3.5 Ice Features 

Sightlines passing through dense molecular clouds typically show additional ab- 
sorption features. The strongest feature is the 3.1 /im O-H stretching mode in 
solid H2O. Whittet et al. (1988) showed that in the Taurus dark cloud complex, 
the strength of this feature is approximately given by 




for Ay 3.3 mag , 

0.093(,4y - 3.3mag) for Ay > 3.3 mag , 



(7) 



Interstellar Dust Grains 



15 



suggesting that dust in regions with Ay > 3.3 has an ice coating, whereas dust 
in regions with Ay < 3.3 is iceless. The 3.1 /im feature is accompanied by a 
number of weaker features due to H2O, plus additional features due to CO2, 
NH3, CO, CH3OH, CH4, and other species. The features due to different species 
often overlap, and the positions and shapes of particular spectral features can 
be sensitive to the composition of the frozen mixture, complicating abundance 
determinations. Table 2 gives the abundances relative to H2O estimated for 
two sightlines into dense star-forming regions (W33A and W3:IRS5) and for the 
sightline to Sgr A* at the Galactic Center. Although all three sightlines have 
similar CO2/H2O ratios, the abundances of other species (e.g., CH3OH) vary 
considerably from one sightline to another. The composition of the ice mantle 
must be sensitive to both local environment and history. 



Table 2: Major Ice Components 



X 


A 


Ar A 




N(X)/N(R 2 0) 




H 


W33A 


W33A 


W3:IRS5 


GC 


H 2 o 


3.1° 


5.5±0.5 d 


1 


1 


1 


XCN 


4.62 


1.4 d 


.035 d 






CO 


4.67 


1.3 e 


.081 e 


.030 e 


< 0.12 h 


H 2 co 


5.8P 


0.17 3 


.065 j 




< .024 h 


CH4 


7.7 


0.11 f 


.016 d 


.0043 / 


.024 h 


NH 3 


9.0 fc 




< 0.05' 


< .035-^ 


0.2-0.3 9 


CH3OH 


9.8 b 


0.9 ± 0.2 d 


0.18 d 


< .004 / 


< .04 h 


C0 2 


15.2 C 


0.58 e 


0.13 e 


0.13 e 


0.14 e 



a Also 6.0, 13.5, 45/tm 

b Also 4.27/tm 

c Also 3.53, 6.85^m 

d Gibb et al. (2000) 

e Gerakines et al. (1999) 

! Giirtler et al. (2002) 

9 Lacy et al. (1998) 

h Chiar et al. (2000) 

j Keane et al. (2001) 

k Also 2.27, 2.96, 3.48/im 

; Taban et al. (2003) 

The grains in star-forming clouds are often aligned (see §6.2), and the radiation 
reaching us from embedded sources is often polarized. Polarization has been 
measured in the features at 3.1 /im (H2O), 4.6 /im (XCN) and 4.67/um (CO) 
(Chrysostomou et al. 1996, Hough et al. 1996). 

3.6 3.3, 6.2, 7.7, 8.6, 1 1. 3 PAH Features 

Many emission nebulae (HII regions, planetary nebulae) and reflection nebulae 
show emission in the 3-15 /im region that is far stronger than expected from 
a > 100A grains heated by the ambient radiation field and emitting thermally. 
Much of the emission is concentrated in five features at 3.3, 6.2, 7.7, 8.6, and 
11.3 /im (see Figure 6). Leger & Puget (1984) identified these features as the 
optically- active vibrational modes of PAH molecules. When H atoms are attached 
to the edge of an aromatic ring skeleton, there are characteristic optically-active 
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Figure 6: PAH emission features in the 5-15 /tm spectrum of the reflection nebula NGC 7023. 
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vibrational modes (Allamandola et al. 1989): 

• C-H stretching mode at 3.3 /im 

• C-C stretching mode at 6.2 /im 

• C-C stretching mode at 7.7/mi 

• C-H in-plane bending mode at 8.6 /jm 

• C-H out-of-plane bending mode with wavelength depending on the number 
of neighboring H atoms: 

— 11.3//m for "mono" H (no adjacent H) 

— 12.0 ^m for "duo" H (2 contiguous H) 

— 12.7//m for "trio" H (3 contiguous H) 

— 13.55 /im for "quartet" H (4 contiguous H) 

Examples of "mono", "duo", "trio", and "quartet" sites are indicated in Figure 
7. Although the broadness of the observed features precludes identification of 
specific PAH molecules, a population of PAH molecules appears to provide a 
natural explanation for the observed emission spectrum, where the vibrational 
excitation is assumed to be the result of "internal conversion" of energy following 
absorption of an optical or UV photon (Leger & Puget 1984, Allamandola et al. 
1989). 

There should in principle be absorption associated with these emission features. 
An absorption feature at 6.2 ^m (see Figure 4) has been observed (Schutte et al. 
1998, Chiar et al. 2000). It matches closely the expected 6.2 //m absorption feature 
(Li & Draine 2001b). 
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Figure 7: The structure of 4 PAH molecules. Examples of mono, duo, trio, and quartet H 
sites are indicated. 

3.7 X-Ray Absorption Edges 

Interstellar grains absorb and scatter X-rays. The X-ray absorption by atoms 
such as C, O, Mg, Si, and Fe will show photoelectric "absorption edges" similar 
to atoms in the gas phase, but the precise energy and structure of the absorption 
edge will be dependent on the chemical nature of the grain material (Forrey et 
al. 1998, Draine 2003b). Significant structure in both scattering and absorption 
is expected near the C K (284 eV), O K (538 eV), Fe L 2 , 3 (708, 721 eV), Mg K 
(1311 eV), Si K (1846 eV), and Fe K (7123 eV) edges (Draine 2003b). The 
Chandra observatory has been used to study the absorption edges of O (Paerels 
et al. 2001, Takei et al. 2002) and O, Mg, Si, and Fe (Schulz et al. 2002) but it 
has not yet proved possible to identify the chemical form in which the solid-phase 
Mg, Si, Fe and O are bound. 

3.8 Extended Red Emission 

The "extended red emission" (ERE) from dust grains provides a potentially im- 
portant clue to the composition of interstellar dust. The ERE consists of a 
broad featureless emission band between ~ 5400A and ~ 9000A, peaking at 
6100A < X p < 8200A, and with a FWHM in the range 600-1000A. The peak 
wavelength X p and band profile appear to vary from one region to another. 

First observed in the Red Rectangle (Schmidt et al. 1980), ERE has since been 
seen in a wide variety of dusty environments, including reflection nebulae (Witt 
& Schild 1985, Witt & Boroson 1990), planetary nebulae (Furton & Witt 1990), 
HII regions (Darbon et al. 2000), the diffuse ISM of our Galaxy (Gordon et al. 
1998, Szomoru & Guhathakurta 1998), and from other galaxies (Pierini et al. 
2002). 
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The ERE must be photoluminescence: absorption of a starlight photon followed 
by emission of a lower energy photoluminescence photon. Based on the detection 
of ERE from the diffuse ISM, the photon conversion efficiency of interstellar dust 
is estimated to be 10%±3% (Gordon et al. 1998, Szomoru & Guhathakurta 1998). 
The photoluminescence efficiency of the ERE carrier must exceed this, since it is 
presumably not the only UV/ visible photon absorber. 

The luminescing substance remains uncertain. Candidate ERE materials in- 
clude HAC (Duley 1985, Witt & Schild 1988), PAHs (d'Hendecourt et al. 1986), 
quenched carbonaceous composite (Sakata et al. 1992), Ceo (Webster 1993), coal 
(Papoular et al. 1996), silicon nanoparticles (Ledoux et al. 1998, Witt et al. 1998, 
Smith & Witt 2002), and carbon nanoparticles (Seahra &: Duley 1999). 

Most candidate materials appear unable to match the observed ERE spec- 
tra and efficiency (Witt et al. 1998). Lab studies of the luminescence spectrum 
of hydrogenated amorphous carbon (HAC) appeared to show a good match to 
observed ERE spectra (Witt & Schild 1988) provided the HAC is annealed suf- 
ficiently to shift the luminescence into the red, where it is observed; however, 
such red-luminescing HAC has a very low photoluminescence efficiency (Furton 
h Witt 1993, Rush et al. 1996) that is incompatible with the high efficiency 
required to explain the ERE from the diffuse ISM. 

Each of the nanoparticle hypotheses appears to have difficulties: (1) Seahra & 
Duley (1999) argue that small carbon particles with mixed sp 2 /sp 3 bonding can 
reproduce the observed ERE profile and required efficiency, but they predict a 
second ERE peak at ~ 1 /j,m which does not appear to be present in NGC 7023 
(Gordon et al. 2000). (2) The PAH hypothesis has difficulty with (a) ERE emis- 
sion and PAH emission having differing spatial distributions in HII regions (Sivan 
&: Perrin 1993, Darbon et al. 2000), and (b) nondetection of ERE in reflection 
nebulae illuminated by stars with T c g < 7000 K (Darbon et al. 1999), whereas 
PAH emission bands are seen in these regions. (3) The silicon nanoparticle hy- 
pothesis appears to be ruled out by nondetection of ~ 20 fim emission that should 
result from single-photon heating (see §9) in the silicon nanoparticle model (Li 
& Draine 2002a). 

Recently, two new emission features have been reported in the near-infrared, 
at 1.15 and 1.5 /im (Gordon et al. 2000). Smith et al. (2001) propose that the 
1.5 /im feature is due to /3-FeSi2. 

Although it would appear that all current proposals have been ruled out, the 
arguments against them are not yet conclusive. For example, perhaps it is possi- 
ble to prepare HAC samples that achieve a higher photoluminescence efficiency, 
or there may be some way to explain the lack of PAH emission seen in the Bub- 
ble nebula, or silicon nanoparticles might be part of larger structures. At this 
point it seems appropriate to continue consideration of HAC, PAHs, and silicon 
nanoparticles as candidate ERE carriers. 

4 PRESOLAR GRAINS IN METEORITES AND INTERPLAN- 
ETARY DUST 

The most primitive meteorites contain presolar grains, which are recognized by 
virtue of isotopically anomalous composition - see Hoppe & Zinner (2000) for a 
recent review. Table 3 lists the six presolar grain types that have been identified 
thus far, along with their abundance (fraction of the total mass) in the bulk 
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meteorite. Nanodiamond (sp 3 -bonded carbon) grains make up 0.05% = 500 ppm 
of the meteorite. A tiny fraction of the nanodiamond particles contain a trapped 
Xe atom; the isotopic abundances of the trapped Xe suggests formation in SN 
ejecta. SiC grains are the next most abundant by mass (6 ppm), with isotopic 
patterns that indicate formation primarily in outflows from AGB stars. About 
0.5% of the SiC grains are of "Type X", with isotopic abundances, including live 
44 Ti (Ti/2 = 59 yr) at the time of grain formation, which suggest formation in 
Type II supernovae. About 1 ppm of the meteorite consists of graphitic (sp 2 - 
bonded carbon) grains, formed in AGB stars, supernovae, and novae. AI2O3 
corundum grains are found, at least some of which appear to have formed in 
outflows from red giants and AGB stars. Finally, there are also S13N4 grains 
which appear to have a SN II origin. 



Table 3: Presolar Grains in Meteorites (Hoppe & Zinner 2000) 



Composition 


diameter ( (im) 


Abundance" 


Origins 6 


C( diamond) 


0.002 


5 x 10" 4 


SN 


SiC c 


0.3-20 


6 x 10~ 6 


AGB 


C(graphite) d 


1-20 


1 x 10 -6 


AGB,SN II,nova 


SiC type X 


1-5 


6 x 10 -8 


SN 


AI2O3 (corundum) 


0.5-3 


3 x 10~ 8 


RG,AGB 


Si 3 N 4 


~1 


2 x 10~ 9 


SN II 



a Overall abundance in primitive carbonaceous chondrite meteorites. 
6 SN = supernova; AGB = asymptotic giant branch star; RG = red giant 
c SiC grains sometimes contain very small TiC inclusions. 
d graphite grains sometimes contain very small TiC, ZrC, and MoC inclusions. 

It is important to recognize that the procedures used to search for presolar 
grains in meteorites begin with chemical treatments designed to dissolve the 
silicate matrix in which the presolar grains are embedded. It therefore comes as 
no surprise that silicates are not present in Table 3. It should also be recognized 
that the formation of even the most primitive meteorites involved temperatures 
high enough to melt rocks, in an environment which was presumably oxidizing 
(since the solar nebula had O/C > 1). Therefore the relative abundances in Table 
3 probably have more to do with the ability of grain types to survive a period of 
exposure to a hot, oxidizing environment than with interstellar abundance. 

Much more informative would be the abundances of presolar grains in cometary 
material. A "core sample" would be ideal, but at this time there is no planned 
mission for returning a sample of pristine cometary ice to Earth. However, the 
Stardust mission (Brownlee et al. 2000) will collect cometary dust on an aerogel 
panel during the 6.1 kms -1 encounter of the Stardust spacecraft with Comet 
Wild 2 in January 2004, and return it to Earth. 

Interstellar grains enter the solar system, and have been detected by the 
Ulysses and Galileo spacecraft; the impact rate is consistent with expectations 
for a k, 0.3 fim grains, but the impact rate for a £ 1 /im grains is much larger 
than expected for the average interstellar dust distribution (Frisch et al. 1999, 
Weingartner & Draine 2001a). The Stardust mission will expose an aerogel panel 
to interstellar grains flowing through the solar system; with ~ 25km s _1 impact 
speeds, these grains are likely to be destroyed, but their residue in the aerogel 
can be analyzed for total mass and composition. We await the return of Stardust 
in January 2006. 
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Active comets release their dust grains, which become part of the interplan- 
etary dust particle (IDP) population. IDPs can be collected from the Earth's 
atmosphere (Brownlee 1985). Bradley (1994) identifies a class of glassy silicate 
IDPs, known as GEMS (Glass with Embedded Metals and Sulfides), which ap- 
pear to be presolar. Most GEMS are between 0.1 and 0.5 fim in diameter - 
similar to the sizes of interstellar silicates (see Figure 8). The 8 — 13 fim infrared 
absorption spectrum of these grains is similar to interstellar silicates (Bradley et 
al. 1999). 

5 CANDIDATE GRAIN MATERIALS 

Based on the observations and clues described above, it is possible to formulate 
a short list of candidate materials for the bulk of interstellar dust. We will omit 
from consideration "rare" materials like AI2O3 - even though Al is depleted from 
the gas phase, given the solar abundance of Al (Al/H = 3.1 x 10 -6 ) aluminum- 
based materials account for only a few % of interstellar grain mass. 

5. 1 Silicates 

There is little doubt that silicate material contributes a substantial fraction of 
the total mass of interstellar dust. As discussed above (§3.2), in the ISM, at least 
~95% of the silicate material is amorphous (Li Sz Draine 2002a). In some cir- 
cumstellar disks the crystalline fraction appears to be higher, although it appears 
that the bulk of the silicates remain amorphous. Bouwman et al. (2001) find that 
5-10% of the olivine mass is crystalline in the dust around Herbig Ae/Be stars. 

Cosmic abundances dictate that the dominant metal ions in silicates will be 
either Mg, Fe, or both. There are two general chemical classes of crystalline 
silicates: Mg-Fe pyroxenes Mg x Fei_ x Si03 (including enstatite MgSi03 and fer- 
rosilite FeSiOs) and olivines Mg2 X Fe2_2a;Si04 (including forsterite Mg2Si04 and 
fayalite Fe2Si04). There are two different crystalline structures for the pyroxenes: 
orthopyroxene and clinopyroxene. 

The optical constants of various crystalline pyroxenes have been measured by 
Jager et al. (1998) and Chihara et al. (2002); crystalline olivines have been mea- 
sured by Jager et al. (1998), Koike et al. (1999), and Suto et al. (2002). 

As noted above, some circumstellar dust exhibits emission features character- 
istic of crystalline silicates, and from the locations and strengths of these features 
it is possible to infer the Mg:Fe ratio. In all cases thus far the crystalline silicate 
material appears to be very Mg-rich and Fe-poor, consistent with pure forsterite 
and enstatite (Tielens et al. 1998, Molster et al. 2002a). Molster et al. (2002b) 
report that enstatite MgSi03 is more abundant than forsterite Mg2Si04 around 
most evolved stars. 

Si is generally heavily depleted in the ISM. If Si is predominantly in silicates, 
and if Fe is not present in silicates, then the interstellar silicates would be expected 
to be enstatite MgSiO-3. The Fe would then presumably be in either metallic Fe 
or in oxide form (FeO, Fe203, FesO^). However, while there is evidence for Mg- 
rich silicate in circumstellar dust, we do not yet have a reliable indication of the 
composition of interstellar silicates, so an appreciable Fe fraction, and olivine-like 
composition, are not ruled out. 

Observed gas-phase abundances of Mg, Fe, and Si, and variations in these 
abundances from one cloud to another, have been interpreted as indicating that 
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silicate grains may have a Mg-rich mantle and Fe-rich core (Spitzer & Fitzpatrick 
1993). The total ratio of (Fe+Mg) atoms to Si atoms in grains is estimated 
to be ~3.4:1 (Fitzpatrick 1997), ~3:1 (Howk et al. 1999) - in excess of the 2:1 
ratio expected for olivines, suggesting that metal oxides or metallic Fe may be 
present in addition to silicates. However, these conclusions depend upon assumed 
values of the total interstellar abundances; Sofia & Meyer (2001) conclude that 
the (Fe+Mg):Si ratio in interstellar dust is actually close to 2:1 for both "halo" 
dust and the dust in the well-studied COph cloud. 

The composition of interstellar silicates remains uncertain, but the study by 
Sofia & Meyer (2001) suggest that olivine MgFeSi04 appears to be a reasonable 
approximation. 

5.2 Carbonaceous Materials 

By "carbonaceous" we mean materials which are predominantly C by mass, in- 
cluding: 

• Diamond: sp 3 -bonded carbon. 

• Graphite: sp 2 -bonded carbon, either monocrystalline or polycrystalline. 

• Amorphous or glassy carbon: a mixture of sp 2 - and sp 3 -bonded carbon, 
with only short-range order. 

• Hydrogenated amorphous carbon: amorphous carbon with an appreciable 
hydrogen fraction. 

• Polycyclic aromatic hydrocarbons: sp 2 -bonded carbon with peripheral H 
(see §3.6). 

• Aliphatic (chainlike) hydrocarbons (see §3.4) 

The presence of sp 2 -bonded carbon is indicated by the 2175A feature and the 
PAH emission bands. The 2175A feature requires sp 2 -bonded carbon in small 
(a ;$ 100A) particles (in larger particles the 2175A feature is suppressed), and 
single-photon excitation of the PAH emission features requires very small par- 
ticles, with ^ 10 4 C atoms. A PAH population incorporating a C abundance 
C/H« 60 x 10~ 6 appears able to account for both the PAH emission features 
(§3.6) and the 2175A feature (§3.1). 

5.3 SiC 

SiC grains are found in meteorites (§4) and an emission feature at 11.3 ^m is ob- 
served in the spectra of many carbon stars (Treffers &: Cohen 1974, Blanco et al. 
1998). However, absence of the 11.3 jiva. feature in the interstellar extinction indi- 
cates that the abundance of Si in SiC dust is less than 5% of the abundance of Si 
in silicates (Whittet et al. 1990). Accordingly, SiC grains, although undoubtedly 
present in the ISM, are not a major component of the interstellar grain mix. 

5.4 Carbonates 

Calcite CaC03 and dolomite CaMg(C03)2 have been detected in dusty disks 
within the planetary nebulae NGC6302 and NGC6537 (Kemper et al. 2002), 
but these carbonates are estimated to contribute less than 1% of the dust mass. 
Carbonates are evidently not a major component of the interstellar grain mix. 
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6 POLARIZATION 

6. 1 Optical- UV Polarization 

Light reaching us from reddened stars is often linearly polarized, because the ex- 
tinction depends on the polarization mode. The degree of polarization as a func- 
tion of wavelength can be approximated by the "Serkowski law" , a 3-parameter 
empirical fitting function (Serkowski 1973): 

p(A) = fmax exp -K(ln(A/A max )) 2 , (8) 

where A max ~ 5500A and K ~ 1 on typical sightlines through diffuse clouds. 
The polarization arises from differential extinction by aligned dust grains, and 
will depend on the degree of alignment with the local magnetic field, the angle 
between the local magnetic field and the line of sight, and the degree to which 
the direction of the magnetic field varies along the line of sight. Studies of many 
sightlines (Serkowski et al. 1975) find that 

Pmax ~ 

0.03A(A max )/mag , (9) 

where A(X) is the (polarization-averaged) extinction at wavelength A. Interstellar 
dust grains must be sufficiently nonspherical and sufficiently aligned so that at 
Amax one polarization mode is extincted ~6% more than the other polarization 
mode. Values of p m axM(A m ax) < 0.03/mag are presumed to result when the 
magnetic field is disordered or not transverse to the line-of-sight, or from regions 
where the degree of grain alignment is lower. 

The wavelength A max varies from one sightline to another, and is correlated 
with the value of Ry, with 

R v » 3.67(A max /5500A) - 0.29 (10) 

(Clayton & Mathis 1988), and 

K » 1.02(A max /5500A) - 0.10 (11) 

( Wilking et al. 1982). The Serkowski law provides quite a good fit to the data for 
0.8 avaT 1 £ A 3 avaT 1 . Observations in the vacuum UV (Anderson et al. 1996, 
Clayton et al. 1996, Wolff et al. 1997) show that the degree of polarization con- 
tinues to decline with decreasing wavelength. This implies that the small grains 
responsible for the rising extinction in the vacuum UV are inefficient polarizers 
- they are either nearly spherical or minimally aligned (Kim &: Martin 1995). 

Most stars show no polarization excess in the 2175A feature but two do: a 
2175A polarization feature is detected in HD 197770 and HD 147933-4 (Wolff et 
al. 1997). 



6.2 Polarized Far-Infrared Emission 

The aligned dust grains that polarize starlight will also produce linearly-polarized 
infrared emission, because for an isotropic dielectric function a submicron grain 
will radiate most effectively when the E field is parallel to the "long" axis of the 
grain. It is now possible to map this polarized emission from dense regions. For 
example, the W51 star-forming region has been mapped at 100 /im (Dotson etal 
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2000), 850 /j,m (Chrysostomou et al. 2002), and 1.3mm (Lai et al. 2001), with 
linear polarizations as large as 10% observed at 1.3mm. 

Hildebrand et al. (1999) have discussed the wavelength-dependence of the mea- 
sured polarization for dense cloud cores and envelopes. The data suggest that 
the grain alignment is larger in warmer regions. 

Observations of quiescent clouds are consistent with no grain alignment in the 
central regions (Goodman et al. 1995, Padoan et al. 2001). 

Radiative torques due to anisotropic starlight play an important role in the 
process of grain alignment (Draine & Weingartner 1997). These torques would be 
ineffective in the central regions of quiescent dark clouds where external starlight 
is heavily attenuated. In star-forming clouds such as M17 or W51, however, 
starlight is provided within the clouds by recently-formed stars. 

7 A PROVISIONAL GRAIN MODEL 

Although we have quite good observational determinations of the extinction in the 
near-infrared, optical, and UV regions, the extinction and scattering properties 
of interstellar grains at very short wavelengths and very long wavelengths are 
probably best obtained by calculating them using a physical grain model that has 
been constrained by a variety of observations, including observations of extinction 
and infrared emission. 

Mathis et al. (1977) discovered that the average interstellar extinction could be 
satisfactorily reproduced by a grain model containing two components - graphite 
grains and silicate grains. Remarkably, the extinction curve was reproduced very 
well if both grain components had power-law size distributions, dn/da oc a~ 3 ' 5 , 
truncated at db 111 1 II 1 IULl 111 S1ZG CL ~ 50A and db 111 8LX 1 11111 111 S1ZG ~ 2500A. The 

success of this grain model led Draine & Lee (1984) to refine the optical constants, 
and extend the treatment into the mid- and far-infrared. The graphite-silicate 
model was applied in the X-ray region by Laor Sz Draine (1993). 

At the time when the "graphite-silicate" model was developed, PAHs had not 
been recognized as a major interstellar grain material. With the recognition 
that the interstellar grain population appears to include a substantial population 
of ultrasmall grains with PAH composition, the "graphite-silicate" model has 
recently been extended in a very natural way: the carbonaceous grains are now 
assumed to be PAH molecules when very small, but to have physical and chemical 
properties that can be approximated by grains of bulk graphite when larger than 
~ 0.01 //m (containing N ^ 10 6 C atoms). The properties of the carbonaceous 
grains are taken to change smoothly from PAH-like to graphite-like as the grain 
size is increased. 

Models developed to account for the infrared emission spectrum of interstel- 
lar dust require that the PAH population contain ~ 15% of the interstellar C 
abundance. These particles necessarily make a major contribution to the UV 
extinction. Weingartner & Draine (2001a) have found size distributions for the 
grains that, including the PAHs, are consistent with the extinction produced by 
the interstellar grain population, and with the size distribution of the PAHs ad- 
justed to reproduce the observed infrared emission from the diffuse ISM (Li & 
Draine 2001b). 

The grain size distribution adopted by Weingartner & Draine (2001a) for aver- 
age Milky Way dust (with Ry = 3.1) is shown in Figure 8, but with abundances 
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reduced by a factor 0.93 to reproduce Af c /Nn from eq. (4). The carbonaceous 
grain distribution is trimodal: The peak at a ~ 0.3 fim is required to reproduce 
the observed extinction curve; the peak at a ~ .0005 fim is required to reproduce 
the 3 — 12 /iin PAH emission features; and the peak at a ~ 0.005 /um improves the 
fit to the observed emission near 60 fim. The peaks at ~ 0.3 /im and ~ .0005 /xm 
are certainly real, but the peak at .005 /um could be an artifact of errors in the 
adopted grain optical and FIR cross sections. 




Figure 8: Size distributions for carbonaceous-silicate grain model of Weingartner & Draine 
(2001a) for Milky Way dust with Rv = 3.1, but with abundances decreased by a factor 0.93 
(see text). 

The dielectric functions of graphite and "astrosilicate" (including structure at 
X-ray absorption edges) have been reestimated by Draine (2003b). The extinction 
and scattering cross sections per H nucleon, calculated for the Weingartner & 
Draine (2001a) grain model, are shown in Figs. 9 and 10. For A ^ 20 /j,m, the 
extinction calculated for this grain model can be approximated by t/Nr ~ 5 x 
10~ 25 cm 2 (A/100/im) -2 , as seen in Figure 10. The scattering varies as expected 
for Rayleigh scattering for A ^ 4/xm: r sca /A^H ~ 1.0 x 10~ 21 cm 2 (A/ //m) . 

This provisional dust model does not include a component to reproduce the 
observed 3.4 //m absorption feature (§3.4). If the graphite in the a }t 200A car- 
bonaceous grains is replaced with a mixture of graphite and aliphatic hydrocar- 
bons, it seems likely that the extinction curve, including the 3.4 fj,m feature, could 
be reproduced with only slight adjustments to the grain size distribution. 

The Weingartner & Draine (2001a) model has Ay/Arg.7 = 14.2 (see Tables 4 - 
6), 25% smaller than the observed value 18.5 for diffuse cloud dust (see Table 1). 
Evidently the silicate abundance or the silicate band strength should be reduced 
by ~ 30%. 
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Figure 9: Extinction and scattering calculated for Weingartner & Draine (2001a) model for 
Rv = 3.1 Milky Way dust, but with abundances reduced by factor 0.93 (see text). 

8 SCATTERING 

8.1 Optical and UV 

Observations of light reflected from dust grains provides another test of dust grain 
models. Individual reflection nebulae are bright, but comparisons with models 
are uncertain because the scattering geometry is in general poorly-determined: 
the illuminating star could be in front of the dust cloud (in which case the typical 
scattering angle would be large) or it could be embedded within it (in which case 
forward scattering by the dust grains between us and the star would dominate). 

The "diffuse galactic light" (DGL) consists of starlight scattered off the dust in 
the diffuse ISM, illuminated by the general interstellar starlight radiation field. 
Since we think we have independent knowledge of the spatial distribution of both 
stars and dust, studies of dust properties using the DGL do not suffer from 
geometric ambiguity. However, the DGL is faint (and must be distinguished 
from the direct light from faint stars), and attempts to infer the grain albedo 
and scattering asymmetry factor g = (cos 6) often find that there can be a one- 
parameter family of (albedo,^) values which are consistent with the observations. 

Figs. 11 and 12 show various observational estimates of albedo and g as a 
function of wavelength for interstellar dust, based on observations of the DGL. 
Determinations of grain scattering properties in the reflection nebulae NGC 7023 
(Witt et al. 1992), IC435 (Calzetti et al. 1995), NGC 2023 (Burgh et al. 2002), 
and the Pleiades (Gibson & Nordsieck 2003) are also shown. 

1. There does appear to be evidence for a decrease in albedo from the optical 
to the vacuum UV, although the results for NGC 7023 and IC435 at A -1 > 
5 yum -1 do not conform. 
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Figure 10: Extinction and scattering calculated for Weingartner & Draine (2001a) model for 
Rv = 3.1 Milky Way dust, but with abundances reduced by factor 0.93 (see text). The dashed 
lines show the asymptotic behavior of the absorption (oc A~ 2 ) and scattering (<x A~ 4 ) cross 
sections. 

2. The data are consistent with (but do not require) a rise in the scattering 
asymmetry factor g from optical to UV. 

3. The data in Figs. 11 and 12 appear to be consistent with the albedo and g 
calculated for the carbonaceous-silicate grain model. 

The scatter in the results of different studies [e.g., the albedo and g near A ~ 
1600A (A" 1 « 6.25 urn' 1 ) in Figs. 11 and 12] serves as a warning that determining 
the albedo is model-dependent, requires accurate photometry, and, in the case 
of photodissociation regions like NGC 7023, NGC 2023, or IC 435, may require 
subtraction of fluorescent emission from H2. 

8.2 X-Ray Scattering 

Interstellar grains scatter X-rays through small scattering angles, and as a result 
X-ray point sources can appear to be surrounded by a diffuse "halo" of X-rays 
scattered by dust grains between us and the source (Overbeck 1965, Hayakawa 
1970, Martin 1970). Measurements of such halos provide a quantitative test of 
interstellar grain models (Catura 1983, Mauche & Gorenstein 1986, Mitsuda et 
al. 1990, Mathis & Lee 1991, Clark et al. 1994, Woo et al. 1994, Mathis et al. 
1995, Predehl k Klose 1996, Smith & Dwek 1998, Witt et al. 2001, Smith et al. 
2002). 

ROSAT observations of Nova Cygni 1992 provide the best observational data 
at this time. Mathis et al. (1995) concluded that the observed halo on day 291 
(after optical maximum) was best explained if the larger (a > 0.1 fim) interstellar 
grains had a "fluffy" morphology, with a void fraction > 25% - if the grains were 
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' Sco 0B2: Gordon et al 1994 [ 
>NCC 7023: Witt et al 1993 
I NGC 7023: Witt etal 1992 
{ G121.0- 10.0: Witt et al 1990 
)L1642: Laureijs et al. 1987 
^NGC7023: Witt etal 1982 
>Coalsack: Mattila 1970 
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4 NGC 2023: Burgh et al. 2002 
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y IC435: Calzetti et al 1995 
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Figure 11: Solid lines: scattering albedo calculated for Weingartner & Draine (2001a) model 
for Milky Way dust with Ry = 3.1, 4.0, and 5.5. Symbols: observational determinations. After 
Draine (2003b). 



compact the scattering halo would be stronger than observed. Smith & Dwek 
(1998) argued that Mathis et al. (1995) had overestimated the scattering due 
to their use of the Rayleigh-Gans approximation, and that in fact the observed 
X-ray halo was consistent with compact grains. The same data was reconsidered 
by Witt et al. (2001), who concluded that the observed halo intensity at small 
scattering angles required a substantial population of large (a > 1 /im) grains. 

Draine & Tan (2003) have revisited this question. Using data for 9 epochs, they 
find that the Weingartner & Draine (2001a) PAH/graphite/silicate grain model 
- with no adjustment to the size distribution - is consistent with the observed 
X-ray halo, for 50"-3000" halo angles. Observations with ROSAT and Chandra 
of other sources at energies between 0.7 and 4 keV also appear to be in general 
agreement with the scattering expected for this grain model (Draine 2003b). 
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Figure 12: Solid lines: Scattering asymmetry factor g = {cos 9) calculated for Weingartner & 
Draine (2001a) model for Milky Way dust with Ry = 3.1, 4.0, and 5.5. Symbols: observational 
determinations. After Draine (2003b). 



9 INFRARED EMISSION 

Heating of interstellar dust grains is primarily by absorption of starlight (colli- 
sional heating dominates only in dense regions in dark clouds - where the starlight 
intensity has been severely attenuated - or in dense, hot, shocked gas). Since the 
heating by starlight photons is quantized and stochastic, the temperature of a 
dust grain is time-dependent. Figure 13 shows the temperature histories of 4 
dust grains, heated by the average starlight background, and cooled by emis- 
sion of infrared photons, over a time span of about 1 day. It is apparent that 
for grains heated by the average starlight background, grains with radii a ^ 200A 
can be approximated as having a steady temperature. Grains with radii a £ 50A, 
however, undergo very large temperature excursions, and the notion of "average 
temperature" is not applicable. Most of the infrared power radiated by such 
small grains occurs during brief intervals following photon arrivals, when the 
grain temperature is close to the peak. 

In order to calculate the time-averaged emission spectrum for such small grains, 
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Table 4: Absorption and Scattering for 5 fim > A > 0.1 fim a 
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a Tabulated data is available at http://www.astro.princeton.edu/~draine/dust/dust.html 
b Cardelli et al. (1989) extinction fit for R v = 3.1 
c Fitzpatrick (1999) extinction fit for R v = 3.1 

d R v = 3.1 Milky Way dust model of Weingartner & Draine (2001a) but with abundances 
reduced by 0.93, and using optical constants from Draine (2003b) 

e A ef f for JHKLL'M from Bessell & Brett (1988); Rjlj from Johnson (1965); Cousins 
UBVRcIc from Fukugita et al. (1996); SDSS ugriz from Gunn (2002, private communica- 
tion). 

f A x = (2.5/lnlO)C7 ext (A)Ar H = 1.086C oxt (A)iVH 

9 ^abs(A) = absorption cross section per unit dust mass. For this model, Mn/M dust = 90. 

one must calculate the energy distribution function dP/dE, for grains of a par- 
ticular size, where dP is the probability that a grain will have vibrational energy 
in interval (E, E + dE). To solve for dP/dE, one must first calculate the specific 
heat for the grain, and the cross sections for absorption and emission of photons. 
Monte-Carlo simulations can be used to obtain dP/dE, but it is more efficient 
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Table 5: Absorption and Scattering" for A > 5 /xm 
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2.36 xlO 2 


25.00 


.0317 


.00031 


-.025 


8.95x10" 


24 


4.79 xlO 2 


18.00 


.0514 


.00066 


-.026 


1.45x10" 


23 


7.76 xlO 2 


14.00 


.0327 


.00283 


-.029 


9.22x10" 


24 


4.92 xlO 2 


12.00 


.0622 


.00301 


-.029 


1.75x10" 


23 


9.37 xlO 2 


9.700 


0.140 


.003 


-.036 


3.94x10" 


23 


2.10 xlO 3 


8.000 


.0491 


.017 


-.050 


1.38x10" 


23 


7.29 xlO 2 


7.000 


.0270 


.055 


-.053 


7.60x10" 


24 


3.85 xlO 2 


6.000 


.0329 


.088 


-.052 


9.27x10" 


24 


4.53 xlO 2 



a See Table 4 

to solve for the steady-state dP/dE (Guhathakurta & Draine 1989, Draine & Li 
2001). 

The notion of "temperature" is somewhat problematic for very small grains, 
but one can define the instantaneous vibrational temperature as being equal to 
the temperature T(E) at which the expectation value of the energy would be 
equal to the actual grain energy; this is the temperature plotted in Figure 13. 
Draine & Li (2001) have shown that a "thermal" estimate for the grain emissivity 
is a good approximation, even for grains containing as few as ~ 30 atoms. 

The infrared and far-infrared glow from interstellar dust can be measured, and 
a dust model must be able to reproduce the observed intensity and spectrum. 
The diffuse ISM at galactic latitudes |6| ^ 20° provides a particularly clean test, 
because we have good estimates for the intensity of the starlight heating the dust, 
and the observed 21 cm emission provides accurate gas column densities which, 
for an assumed dust-to-gas ratio, determine the column density of dust. Although 
the surface brightness of the infrared emission is low, it has been detected by IRAS 
and CO BE after averaging over large areas of sky. Figure 14 shows the observed 
emissivity per H nucleon, based on COBE-FIRAS and COBE-DIRBE photome- 
try. 3.0-4.5 /im and 4.5-11.3 /im spectra, measured at b ~ by the IRTS satellite 
have been added to the plot, normalized to agree with the DIRBE photometry 
in the 3.5, 5.0, and 12 fim bands. 

In this dust model, approximately 20% of the absorbed starlight energy is 
reradiated at A < 12 ^m by PAHs. The same grain model, applied to reflection 
nebulae, can acount for the observed dependence of dust emission spectrum on 
temperature of illuminating star (Li & Draine 2002b). With appropriate changes 
to the grain mixture, the model can also account for the observed IR emission 
from the SMC (Li & Draine 2002c). 

Approximately two thirds of the radiated power is at A ^ 50 /im - this is emis- 
sion from a £ .01 /ira grains that are maintained at a nearly steady temperature 
~ 15 — 20 K by starlight. The remaining ~ one third of the radiated power is 
from grains with radii a ^ 50A that are cooling following "temperature spikes" 
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Table 6: Absorption and Scattering" for A < 0.1 um 



A 


his 




albedo 


(cos(0)) 


s~i a 
ext 




K a 
abs 




eV 


model 






cm 2 /H 




cm 2 /sf 


.09120 


13.6 


8.256 


0.248 


0.658 


2.33 x 10" 


21 


9.37x 10 4 


.04881 


25.4 


5.574 


0.291 


0.813 


1.57x10" 


21 


5.97x 10 4 


.02339 


53.0 


3.355 


0.366 


0.953 


9.46 x 10" 


22 


3.21 x 10 4 


.02279 


54.4 


3.577 


0.369 


0.949 


1.01x10" 


21 


3.41 xlO 4 


.02214 


56.0 


3.469 


0.393 


0.950 


9.78 x 10" 


22 


3.18x 10 4 


.01240 


100 


2.564 


0.452 


0.988 


7.23 x 10" 


22 


2.12x 10 4 


.01170 


106 


2.908 


0.432 


0.985 


8.20 x 10" 


22 


2.49 x 10 4 


.006199 


200 


2.000 


0.578 


0.997 


5.64 x 10" 


22 


1.27x 10 4 


.004397 


282 


1.426 


0.619 


0.999 


4.02 x 10" 


22 


8.20 x 10 3 


.004350 


285 


2.355 


0.452 


0.998 


6.64 x 10" 


22 


1.95 x 10 4 


.004133 


300 


1.812 


0.517 


0.999 


5.11 x 10" 


22 


1.32 x 10 4 


.002480 


500 


0.971 


0.632 


1.000 


2.74 x 10" 


22 


5.40 x 10 3 


.002317 


535 


0.691 


0.386 


1.000 


1.95x 10" 


22 


6.40 x 10 3 


.002305 


538 


1.325 


0.355 


1.000 


3.74 x 10" 


22 


1.29 x 10 4 


009996 


540 


1.138 


0.453 


1.000 


3 91 x 1 0" 


22 


9 40 x 1 3 


.002254 


550 


1.056 


0.478 


1.000 


2.98x10" 


22 


8.32 xlO 3 


.001759 


705 


0.575 


0.433 


1.000 


1.62x10" 


22 


4.92 xlO 3 


.001751 


708 


1.219 


0.353 


1.000 


3.44x10" 


22 


1.19xl0 4 


.001734 


715 


0.799 


0.515 


1.000 


2.25x10" 


22 


5.84 xlO 3 


.001720 


721 


0.941 


0.432 


1.000 


2.65x10" 


22 


8. 07 xlO 3 


.001710 


725 


0.845 


0.503 


1.000 


2.38x10" 


22 


6.33 xlO 3 


.001550 


800 


0.751 


0.514 


1.000 


2.12x10" 


22 


5.51 xlO 3 


.001240 


1000 


0.563 


0.564 


1.000 


1.59x10" 


22 


3.71 xlO 3 


.000620 


2000 


0.181 


0.584 


1.000 


5.10x10" 


23 


1.14 xlO 3 



a see Table 4 



resulting from absorption of individual starlight photons. We can immediately 
see the following: (^4) The grain population must contain enough a ^ 50A grains 
to account for ~ one third of the absorption of starlight. (B) The resemblance of 
the 3-12 fim emission spectrum from the diffuse ISM to the emission spectrum 
expected from polycyclic aromatic hydrocarbon (PAH) molecules, radicals, and 
ions implies that a large fraction of the smallest grains must be PAHs. Figure 
15 shows the emission calculated for the PAH/carbon/silicate dust model illumi- 
nated by the average interstellar radiation field spectrum. 

The far-infrared emission is a good tracer of interstellar dust. The IRAS and 
DIRBE all-sky maps can be used to infer the dust column density (Schlegel et 
al. 1998). 

10 MICROWAVE EMISSION 

10.1 Spinning Grains 

Theoretical models to account for the observed infrared emission require large 
numbers of extremely small dust grains that, following absorption of starlight 
photons, radiate in their optically-active vibrational modes. In addition to vi- 
brating, these grains will rotate, with rotation rates reaching tens of GHz. If they 
have electric dipole moments, there will be rotational emission from these grains 
at microwave frequencies. Various workers have considered microwave emission 
from spinning grains (Erickson 1957, Hoyle k. Wickramasinghe 1970, Ferrara h 
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Figure 13: A day in the life of 4 carbonaceous grains, heated by the local interstellar radiation 
field. r a b s is the mean time between photon absorptions (Draine & Li 2001). 



Dettmar 1994). The rotational emission from this grain population has been es- 
timated by Draine & Lazarian (1998a, b), who showed that it could account for 
dust-correlated microwave emission that was discovered by Kogut et al. (1996a,b) 
in studies of the angular structure of the cosmic microwave background. 

The theoretical estimate for the rotational emission, plus the low-frequency 
tail of the vibrational emission, is shown in Figure 16. Figure 16 also shows 
results from a number of observational studies - some over large areas of the 
sky, others over smaller regions. Finkbeiner et al. (2002) detect dust-correlated 
microwave emission in pointed observations of a dark cloud (L1622) and toward 
LPH 201.663+1.643, a diffuse HII region behind heavy obscuration. Most re- 
cently, Finkbeiner et al. (2003a) use 8 and 14 GHz surveys to map the anomalous 
microwave emission from the Galactic plane. 

With the exception of the observations toward LPH 201.663+1.643, the ob- 
served microwave emission appears to be generally consistent with the theoretical 
prediction for emission from spinning dust grains. 

The observed microwave emission from LPH 201.663+1.643 is, however, an 
order of magnitude stronger than would be expected from the observed column 
of dust. While the rotational emission from spinning dust grains does depend 
on environment Draine & Lazarian (1998b), it does not appear possible to ex- 
plain the strong emission observed in this direction using reasonable variations 
in environmental conditions. 



10.2 Magnetic Grains? 

Although rotational emission from spinning grains is a natural prediction for a 
dust model with large numbers of small grains, other emission mechanisms are 
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- Observed Emission from Interstellar Dust - 




Figure 14: Observed emission from diffuse interstellar dust. Crosses: IRAS (Boulanger & 
Perault 1988); Squares: COBE-FIRAS (Finkbeiner et al. 1999); Diamonds: COBE-DIRBE 
(Arendt et al. 1998); Heavy Curve: IRTS (Onaka et al. 1996, Tanaka et al. 1996). The total 
power ~ 5.1 x 10~ 24 ergs _1 /H is estimated from the interpolated broken line. 



also possible, including thermal magnetic dipole radiation from magnetic grains 
(Draine &: Lazarian 1999). We know that Fe contributes a significant fraction of 
the interstellar dust mass; if some of this Fe is in ferrimagnetic or ferromagnetic 
compounds, grains containing these compounds would emit magnetic dipole ra- 
diation. The idea is that the magnetization of the magnetic domains in the grain 
would fluctuate around the minimum-energy state; these fluctuations would be 
at GHz frequencies, and the fluctuating magnetization would generate electro- 
magnetic radiation. If an appreciable fraction of the Fe in interstellar grains were 
in magnetite Fe30"4, it is estimated that the resulting thermal magnetic dipole 
emission would be comparable to the observed microwave emission. The power 
radiated per grain is simply proportional to the mass of magnetic material in that 
grain, and proportional to the grain temperature. Since it is presumed that most 
of the magnetic material would be in large grains, it would be expected that the 
large grains would account for the bulk of the magnetic dipole emission. 

It is therefore possible that the observed microwave emission may contain ap- 
preciable contributions from both (ultrasmall) spinning grains and (large) mag- 
netic grains. 

11 ABUNDANCE ISSUES 

Dust population models that are constructed to reproduce the observed wavelength- 
dependent extinction are found to imply substantial abundances of elements in 
grain material - approaching or exceeding the abundances believed to be appro- 
priate to interstellar matter. Kim & Martin (1996) referred to the "C/H crisis", 
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Figure 15: Solid line: infrared emission calculated for WD01 Rv = 3.1 Milky Way dust model 
(but with abundances reduced by factor 0.93) for grains heated by diffuse starlight (Li & Draine 
2001b). Squares: COBE-FIRAS (Finkbeiner et al. 1999) Diamonds: COBE-DIRBE (Arendt et 
al. 1998). 

and Mathis (1996) has argued that the problem may be even more severe for Mg, 
Si, and Fe required for silicate grains. 

Whether there is a problem or not depends on the adopted values for to- 
tal interstellar abundances. Mathis (1996) took the position that interstellar 
abundances should be taken to be ~70% of solar. However, recent changes in 
inferred solar abundances for C and O, and comparison of photospheric abun- 
dances in F and G stars with those in B stars, led Sofia & Meyer (2001) to reject 
the hypothesis that interstellar abundances are appreciably subsolar. Sofia &: 
Meyer estimate the dust-phase abundances on the sightline to £Oph if interstel- 
lar abundances are taken to be the same as photospheric abundances in F and 
G stars. From these abundances, we can construct carbon grains with volume 
= 2.05 x 10 _27 cm 3 /H (for an assumed p = 2.2 gem" 3 ), and silicate/Si02 grains 
with volume 2.90 x 10 _27 cm 3 /H (for an assumed p = 3.5 gem -3 ). 

The provisional grain model of §7 - the WD01 grain model with abundances 
reduced by a factor 0.93 - requires a carbonaceous grain volume V^ ar b = 2.10 x 
10~ 27 cm 3 /H, and a silicate grain volume V s n = 3.68 x 10~ 27 cm 3 /H: thus we see 
that the carbon abundance is in line with observations, but the provisional grain 
model requires ~26% more Mg, Si, and Fe than appears to be available. 

Is this 26% discrepancy a serious problem? Consider the following: 

1. 6 of the 14 sightlines in Figure 2 have Aj/Nr differing by more than 30% 
from the nominal value. This may represent measurement errors, but sug- 
gests real variations in Aj/N^. Since it is not clear how the sample of 
sightlines is constructed, it is not apparent that the average value of Nn/Aj 
for the sample is the same as the overall dust-to-gas ratio for the ISM. 
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Figure 16: Dust-correlated microwave emission, expressed as microwave surface brightness per 
unit column density of hydrogen. Curves: theoretical estimate of Draine & Lazarian (1998a,b). 
Symbols: observational determinations of "anomalous microwave emission". 



2. The dust model assumes the grains to be spheres. Real grains are nonspher- 
ical. This will tend to provide somewhat more extinction per unit mass of 
grain material; this could increase the visual extinction per unit grain mass 
by ~ 10 - 20%. 

3. The dust model assumes the grains to be compact, but real grains may 
contain voids. This could conceivably increase the visual extinction per 
unit grain mass by ~ 10%. 

4. The dust model uses "bulk" optical constants and Mie theory. The outer 
few monolayers of grain material, however, will have electronic structure 
(and therefore optical properties) differing from "bulk" material, so that 
the model calculations should not be considered to be precise even if the 
grains were actually spheres and we had the correct bulk optical constants. 

At the present stage of observational and theoretical knowledge, we do not con- 
sider the apparent "abundance crisis" to be sufficient grounds to reject our pro- 
visional grain model, but new data could change this conclusion. 

12 GRAIN ALIGNMENT 

The polarization of starlight (§6) was discovered over 50 years ago (Hall 1949, 
Hiltner 1949), and was immediately recognized as being due to aligned dust 
grains. The mechanism responsible for the alignment has been sought by theorists 
for the past half century. Two separate alignments are involved: (1) alignment 
of the grain's principal axis of largest moment of inertia with the grain's angular 
momentum J, and (2) alignment of J with the galactic magnetic field. 
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The alignment processes turn out to involve quite subtle physics [for an intro- 
duction, see Draine (2003a)]. The paramagnetic dissipation process proposed by 
Davis & Greenstein (1951) would be able to gradually align grains, but recent 
work shows that radiative torques due to starlight can be much stronger than the 
Davis-Greenstein torque; the starlight torques can affect both the rate of grain 
rotation - driving grains to suprathermal rotation rates - as well as the direction 
of the grain angular momentum (Draine & Weingartner 1996, 1997). The elec- 
tronic and nuclear paramagnetism of the grain couple the rotational degrees of 
freedom to the vibrational degrees of freedom, with important consequences for 
the rotational dynamics of a tumbling grain (Lazarian & Draine 1999a, b). 

H2 formation on grain surfaces is thought to take place at a limited number 
of reaction sites, and the recoil from nascent H2 provides a systematic torque on 
the grain that is fixed in body coordinates (Purcell 1979), and that can drive a 
grain to suprathermal rotation rates. When the grain rotational kinetic energy 
is not too large compared to kT^, where is the grain temperature, the process 
of "thermal flipping" (Lazarian & Draine 1999a) can cause the H2 formation 
to time-average to zero, so that suprathermal rotation does not take place - a 
process termed "thermal trapping". Radiative torques due to starlight, however, 
are not fixed in grain body coordinates and therefore do not time-average to zero 
even when thermal flipping is rapid, so that grains subject to starlight torques are 
immune to "thermal trapping" , and are expected to rotate suprathermally most 
of the time. Starlight torques are important only for grains with radii a ^ 0.1 /zm. 
This may explain the fact that grains with radii a ^ 0.2 //m appear to be much 
more strongly aligned than grains with radii a £ 0.05 fim (Kim & Martin 1995). 

The theory of grain alignment is complicated by the need to describe the dy- 
namics of tumbling grains during intervals when the grain rotational kinetic en- 
ergy is of order kT gas , including the effects of thermal flipping (Weingartner & 
Draine 2003). Progress is being made, and quantitative calculations of grain 
alignment for model grains should be possible within the next year or two - 
this should allow us to "predict" the wavelength-dependent polarization for a 
grain model once a geometric asymmetry is assumed for the grains. At that time 
we will learn whether our current picture of grain dynamics is consistent with 
observations. 

The very smallest grains spin with GHz rotation rates and emit electric dipole 
radiation in the microwave (§10). If these grains are aligned with the magnetic 
field, the rotational emission will be polarized. It appears, however, that these 
grains will be only minimally aligned, with degrees of alignment of only a few 
percent (Lazarian & Draine 2000). The Davis-Greenstein paramagnetic dissipa- 
tion mechanism is suppressed in these grains by the unavailability of vibrational 
modes with energies TiuJ TO t- 

13 EVOLUTION OF THE GRAIN POPULATION 

The interstellar grain population in a galaxy reflects the interplay of many pro- 
cesses. 

Are interstellar grains "stardust"? There is ample evidence that grains form 
in some stellar outflows - we see infrared emission from dust in outflows from 
red giants and carbon stars, and in planetary nebulae. In all of these cases it is 
believed that the bulk of the dust condensed out of initially dust-free gas from the 
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stellar envelope. The 10 /jm silicate feature is seen in outflows from oxygen-rich 
stars, consistent with the expectation that silicates should be a major condensate 
in cooling gas with O/C > 1. The 10 //m feature is not seen in outflows from 
stars with C/O > 1, where instead one sometimes sees a feature at 11.3 /im 
that is attributed to SiC (Treffers & Cohen 1974, Blanco et al. 1998), consistent 
with the expected condensates from cooling gas with O/C < 1. Because isotopic 
anomalies are used to find them, presolar grains in meteorites generally appear 
to be "stardust" - a given grain was formed in the outflow/ejecta from a single 
star with a particular isotopic composition. 

There is no doubt that dust is formed in many stellar outflows. However, this 
does not imply that the bulk of the dust in the ISM was produced in this way. 

The problem is that interstellar grains are subject to sputtering (Draine & 
Salpeter 1979a) and grain-grain collisions in interstellar shock waves, and the- 
oretical studies of grain destruction in the Milky Way find a mean "residence 
time" of an atom in a grain of order Td es t ~ 3 x 10 8 yr (Draine & Salpeter 1979b, 
Jones et al. 1994). 

Now suppose that all Si atoms enter the ISM as constituents of dust grains. 
The total mass of the ISM is Mjsm ~ 5 x 1O 9 M0, and star formation takes place 
at a rate M* < 5M Q yr" 1 . Thus the mean residence time of a metal atom in the 
ISM is tism = Mism/M* ss 10 9 yr. Only a fraction r dcst /(r dcst + tism) ~ 0.2 of 
the Si atoms in the ISM would be found in the original stardust particle in which 
they left a star. 

However, when we observe gas-phase abundances in the ISM, we typically find 
£ 90% of the Si atoms missing from the gas phase. Unless we have grossly 
underestimated the grain destruction time Td e st> it must follow that most of the 
Si atoms in interstellar grains joined the grain in the ISM! 

Therefore most interstellar dust is not stardust - the composition of interstellar 
grains is determined by physical processes in the ISM. The importance of inter- 
stellar processing is also evident from the large variations seen in the extinction 
curve from one sightline to another - these variations must be the result of in- 
terstellar effects, since any given sightline would be expected to average over the 
stardust produced by many individual stars. 

The interstellar processes that alter grains include 

• grain erosion by sputtering; 

• grain shattering in low-velocity grain-grain collisions; 

• grain vaporization in high-velocity grain-grain collisions; 

• grain coagulation; 

• accretion of impinging atoms and molecules on grain surfaces; 

• chemical reaction of impinging atoms and radicals with grain surfaces; 

• photolysis of grain material by UV starlight; 

• alteration of grain material by cosmic ray/X-ray irradiation; 

• photodesorption of atoms and molecules from grain surfaces. 

In addition, the local grain mixture can be modified by size- and composition- 
dependent drift of dust through the gas, driven by anisotropic starlight: relative 
drift speeds ^ 0.1km s^ 1 can displace large grains from small grains by ~few pc 
in 10 7 yr (Weingartner & Draine 2001b). 
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Understanding the balance among these different processes requires not only 
an understanding of the microscopic physics but also a quantitative model for 
the transport of grains from one "phase" of the ISM to another. Because grain 
destruction appears to be dominated by processes associated with interstellar 
shock waves, the propagation of these shock waves through the inhomogeneous 
ISM must be modeled. 

The required modeling is complex, and we are not yet at a stage where "a 
priori" modeling of the evolution of the interstellar grain population can be con- 
sidered reliable. However, the observed depletions can be used to constrain mod- 
els of the dynamic ISM: in order for grains to be able to "scavenge" elements 
like Si from the gas, there must be rapid "cycling" of matter between the diffuse 
low-density phases (where most grain destruction takes place) and denser regions 
(where densities are high enough for grain growth to be effective), with cycling 
times as short as ~ 10 7 yr (Draine 1990). 

14 FUTURE DIRECTIONS 

Our present understanding of interstellar grains is built upon an observational 
base. New observational facilities can be expected to produce important exten- 
sions to this observational base, possibly requiring revisions to our models. 

• SIRTF (Space InfraRed Telescope Facility) will provide low and medium 
resolution 5-40 /mi spectra of reflection nebulae, 55-96 /im spectral energy 
distributions, and imaging in 7 photometric bands from 3.5 to 160 /mi. We 
can anticipate detection of new emission features in the 15-40 /mi region 
(where ISO had limited sensitivity), and studies of variations of dust emis- 
sion features and spectra with environment. 

• Astro-F (Pearson et al. 2002 and references therein), scheduled for launch 
in 2004, will obtain all-sky maps in 4 bands between 50 and 200 /im, with 
improved spatial resolution and much better sensitivity than IRAS. This 
will greatly improve our knowledge of the distribution of diffuse dust, the 
degree to which dust is correlated with H I 21cm emission, and spatial 
variations in the far-infrared band ratios. Astro-F may be able to provide 
all-sky imaging in 6-11 /mi that would provide a map of PAH emission. 
Correlation of the PAH emission with microwave emission would test the 
hypothesis that the microwave emission comes from spinning grains. 

• WMAP (Wilkinson Microwave Anisotropy Probe) is providing maps of the 
microwave sky in 5 bands between 22 and 90 GHz, with FWHM of 0.53° at 
40 GHz. These maps will help characterize the spectrum of the microwave 
emission from dust, as well as variations in dust emission (intensity and 
spectrum) from one region to another (Bennett et al. 2003). It will be of 
great interest to cross-correlate the all-sky microwave maps with infrared 
and far-infrared maps from IRAS, COBE, and, when available, Astro-F. 

• Planck/Surveyor, scheduled for launch in 2007, will provide all-sky submm 
maps that will determine the diffuse dust emission spectrum from 30 - 857 
GHz (1 cm - 350 /im). 

• Herschel, scheduled for launch together with Planck, will be capable of 
sensitive broadband diffraction-limited imaging at 60-700 /im, extending 
beyond the 160 /mi cutoff of SIRTF. The longer wavelength coverage will 
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provide improved determinations of the emission from cool dust, and will 
allow a census of dust masses in nearby galaxies. 

• The Stardust spacecraft is expected to return to Earth in 2006. Its aerogel 
panels will have interstellar grain material at impact sites and (we hope) 
intact cometary dust particles. 

• SOFIA (Stratospheric Observatory for Infrared Astronomy) will provide 
higher angular resolution dust emission imaging between 5 and 240 fim 

• Ground-based millimeter and submillimeter telescopes and arrays (e.g, BIMA, 
JCMT, CSO, ALMA) can map thermal dust emission from dense regions, 
including polarimetric mapping. 

Many problems await solution, some of which are highlighted here: 

• The carrier of the 2175A feature remains unidentified. The PAHs required 
by the observed IR emission features are expected to produce strong absorp- 
tion in this region, but this remains conjecture. We require lab measure- 
ments of UV absorption by PAH molecules, radicals, and ions containing 
20-10 3 C atoms for comparison with the observed 2175A profile. 

• The nature of the ERE carrier remains an open question. The high in- 
tensities reported for diffuse clouds imply that the ERE carrier must be a 
major grain constituent (§3.8). Further observations of the ERE are needed, 
particularly to confirm the intensity of ERE from diffuse clouds. Further 
spectroscopy of ERE emission from bright regions, and laboratory studies 
of candidate ERE carriers, may lead to identification of the ERE carrier. 

• We should be able to identify at least some of the DIBs! Gas-phase spec- 
troscopy of candidate molecules, radicals, and ions is needed. When the 
right candidates are studied, positive spectroscopic identification of DIBs 
will follow. 

• PAH molecules, radicals, and ions appear to be a major part of the grain 
population in the Milky Way, but the origin of these PAHs is unknown. 
Further observational study of spatial variations in PAH abundance (and 
perhaps spectral properties) may provide valuable clues. 

• We would like to understand the balance of processes (growth, coagulation, 
erosion, shattering) accounting for the size distribution of interstellar dust. 
Observational determinations of regional variations in the size distribution 
will frame the theoretical attack on this problem. 

• Are there really separate populations of carbonaceous grains and silicate 
grains? If so, how does grain growth in the ISM maintain these separate 
populations? 

• Is the anomalous microwave emission due primarily to spinning dust grains? 
Is there a significant contribution from magnetic dust materials (§10)? 

• Can we diagnose the composition of interstellar dust using X-ray absorption 
measurements(§3.7)? 

• The problem of interstellar grain alignment (§12) seems ripe for solution. 
Radiative torques due to starlight appear to be an important part of the 
alignment process. A solution to this problem may be in sight. 
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